Generalized Berreman's model of the elastic surface free energy of a nematic liquid 

crystal on a sawtoothed substrate 
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In this paper we present a generaUzation of Berreman's model for the elastic contribution to the 
surface free-energy density of a nematic liquid crystal in presence of a sawtooth substrate which 
favours homeotropic anchoring, as a function of the wavenumber of the surface structure q, the 
tilt angle a and the surface anchoring strength w. In addition to the previously reported non- 
analytic contribution proportional to —qlnq, due to the nucleation of disclination lines at the wedge 
bottoms and apexes of the substrate, the next-to-leading contribution is proportional to q for a 
given substrate roughness, in agreement with Berreman's predictions. We characterise this term, 
finding that it has two contributions: the deviations of the nematic director field with respect to a 
reference field corresponding to the isolated disclination lines, and their associated core free energies. 
Comparison with the results obtained from the Landau-de Gennes model shows that our model is 
quite accurate in the limit wL > 1, when strong anchoring conditions are effectively achieved. 

PACS numbers: 61.30.-v,61.30.Dk,61.30.Hn,61.30. Jf 
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I. INTRODUCTION 

The behaviour of nematic hquid crystals in the pres- 
ence of microstructured substrates has been subject of 
intensive research in the recent times [H-il. This problem 
has practical apphcations such as the design of zcnithally 
bistable devices or the trapping of colloidal parti- 

cles in specified sites [1, [l^ . It is well known that the 
nematic director field, in presence of the structured sub- 
strate, may be distorted, leading to an elastic contribu- 
tion to the free energy. Since the seminal work by Berre- 
man [ll|, [T^], this problem has been extensively stud- 
ied and generalized in the literature 0, IT3i - [26| . Wetting 
and filling transitions by nematic on these grooved sur- 
faces have also been studied [27l - l30| . When the substrate 
presents cusps, topological defect nucleate on them (Tsl - 
[isl . [25I . [26j , and Berreman's expression of the elastic con- 
tribution to the free-energy density, which is proportional 
to the wavenumber of the substrate structure q, breaks 
down because of the emergence of a non-analytical con- 
tribution pro por tional to —qlnq associated to the nucle- 
ated defects [2^. This result constrasts with the phe- 
nomenology observed in smooth substrates, as sinuosidal 
substrates, where in absence of topological defects a suit- 
able generalization of Berreman's model works [23l. [soj. 

In this paper we will complete the characterization of 
the elastic contribution to the surface free-energy den- 
sity for sawtoothed substrates [2^. Beyond the —qlnq 
term previously mentioned, we find that the next-to- 
leading contribution follows Berreman's scaling with the 
wavenumber q. The origin of this term is twofold: the 
deviations of the nematic director field with respect to 
the distortions imposed by the presence of the nucleated 
disclination lines, and the defect core contributions. We 
estimate both contributions, finding a fairly good agree- 
ment with the reported values in Ref . [2^ . 

The paper is organized as follows. The problem is set 



up in Section II, where we identify the different contribu- 
tions to the elastic contribution to the surface free-energy 
density. Sections III and IV are devoted to the estima- 
tion of these contributions, and the obtained results will 
be discussed in Section IV. We will end up with the con- 
clusions in Section V. 



II. THE MODEL 



We consider a nematic liquid in contact with a saw- 
tooth substrate characterized by the angle a and the 
length side L (see Fig. [T]). The substrate favours 
homeotropic anchoring of the molecules. We assume 
traslational symmetry along z axis and a periodic distri- 
bution of wedges and cusps along the x axis. Under these 
conditions, the nematic director field n(r) shows only in- 
plane distortions , so it can be parametrized by using 
the angle between the local director and the y axis, 
yielding n(r) = (— sin6'(r), cos 6'(r), 0). Far from the 
substrate, the bulk nematic phase orients homogeneously 
along either the x axis (perpendicular texture N-^) or the 
y axis (parallel texture TV"). The nematic order of the 
system can be locally represented by a traceless symmet- 
ric order parameter second-rank tensor Q, with Cartesian 
components Qij = ^S[n.inj — ^Sij] + ^B[lilj — mimj], 
where S is the nematic order parameter, which measures 
the orientational ordering along the nematic director, and 
B the biaxiality parameter, which measures the ordering 
of the molecules on the orientations perpendicular to n, 
characterized by the eigenvectors 1 and m. 



As in previous works |26|, |28|, |29| , the system will be de- 
scribed within the Landau-de Gennes (LdG) framework. 
The order parameter tensor profile is obtained by mini- 
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mizing the surface free-energy density functional 



At coexistence r = 1, and Eq. ([2]) reduces to: 
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where A = 2L cos a is the substrate periodicity wave- 
length, Q = Q/Sb, where Sb = Sb{T) is the bulk value 
of the nematic order parameter at the temperature T, 
and T = Sb{T) / Sb{TN i) is the ratio between the ne- 
matic order parameter at the temperature T and the 
nematic- isotropic transition temperature Tni- The po- 
sitions are measured in units of v^^, where ^ is the 
nematic correlation length along the local nematic di- 
rector axis. Finally k is the ratio between the rele- 
vant elastic parameters (k > —3/2) and w the (reduced) 
anchoring strength. Homeotropic alignment of the ne- 
matic is favoured by setting Qsurf = (Si' CS) — l)/2, 
being v = ((x/lxl) sin a, — cos a, 0) the outwards nor- 
mal vector to the substrate. The global minimum of 
the functional Eq. ([T|) yields the mean-field equilib- 
rium surface free-energy density, /, in appropiated units. 
The contribution due to the elastic deformations induced 
by the substrate structure, /elastic, can be obtained as 
/clastic = / — raNwiw), where the surface roughness is 
r = 1/ cos a and aNW (w) is the nematic- flat substrate 
interfacial tension. In a similar way as it was obtained at 
nematic isotropic coexistence in Ref. (soj . the interfacial 
tension for the LdG model has the expression: 



where ^(0) = (1 + \/l + 2V2w)/2. 

A systematic study of this system via full minimiza- 
tion of the LdG functional was done in Ref. [gg]. For 
this purporse, the functional Eq. ^ was numerically 
minimized by using a conjugate-gradient method. The 
numerical discretization of the continuum problem is per- 
formed with a finite element method [3]| combined with 
adaptive meshing in order to resolve the different length 
scales that may emerge in the problem [s^l • It was found 
that the N-^ texture has lower free energy if a < 7r/4 
owing to lesser distortion. Conversely, the A^" texture 
has lower free energy for a > 7r/4, in agreement with 
earlier predictions [isl - fTEi [25|. For large wL, strong an- 
choring conditions are effectively achieved, leading to the 
nucleation of disclination lines characterized by non-half- 
integer winding numbers along the ridges and wedges of 
the substrate [isl - fTsl . [25l . [26j . As a consequence, the elas- 
tic contribution to the surface free-energy density has the 
following scaling [2^ 

/C(a) 9 cos a 
/clastic « 77— <1^^ — Z ^~7^^(", w;) (5) 



2tt 



TT 2TT 



where q — 2t:/X ~ n/L cos a is the substrate periodicity 
wavenumber, IC{a) is defined as: 

, -^J"' texture 
(f ) 

/C(a) = { ^ (6) 
Knt^ A^ll texture 



being K = (9/2)(2-f k)/(3 + 2k) the reduced bulk elastic 
constant associated to bend and splay distortions. From 
the numerical results, the function B is found to depend 

cnw = ' ;r^^ — , , on the substrate rougness (i.e. a) and nematic texture, as 

well as the anchoring w, but asymptotically not on L for 
large wL. However, as L increases, the complete min- 
imization becomes very time-consuming. On the other 
hand, we do not get information about the origin of 
B{a, w). 

In this paper we will introduce an alternative way 
to obtain /elastic from the functional Eq. ([1]). We di- 
vide the minimization domain into three regions (see Fig. 
t5'(0) + T-1 + \/7^(r - 1)(S'(0) + 1)2 + 5(0)2 \[I]): most of the domain wih correspond to the region /, 
X In I ■_ _ . = I formed by the points which are far enough from the sub- 

strate. The neighbourhood of the substrate will be split 
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where 5(0) = 5(0)/5f, and 5(0) is the nematic order 
parameter at the substrate. The value of 5(0) is obtained 
as the largest solution of the equation: 

(3-2T)52(0)-253(0) + r5'*(0) = (3) 



into two regions: the region //, formed by the union 
of the circular sections of radii £^ < Rc L centered 
at each wedge and apex; and region ///, which are the 
points which are at a distance smaller than 77 ^ ^ to 
the substrate, but at distances larger than Rc from any 
substrate ridge. Our hypothesis is that, for large wL, 
the minimization of the surface free-energy functional re- 
stricted to each region (subject to appropriated boundary 
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conditions), gives a good approximation to the complete 
minimization of T . On the other hand, we anticipate 
that this analysis will give us some insight in the differ- 
ent contributions to B{a,w). 

We start with the minimization of region ///. As it 
was argued in Ref. [l^, large wL leads to strong an- 
choring conditions. So, in order to minimize the surface 
free-energy density, we impose to the angle field 9{r) to be 
constant along its boundary, so the nematic director field 
is homogeneous and equal to the normal to the substrate. 
Consequently, the minimization of the free-energy func- 
tional in this region will lead to a homogeneous director 
field normal to the substrate, although the nematic order 
parameter S at each point will depend on its distance to 
the substrate. We impose the following fixed boundary 
conditions for S: S = St at the boundary between regions 
I and and the equilibrium nematic order parameter 
profile for the flat wall case at the boundary between 
regions // and ///. Assuming that rj is large enough, 
this situation is completely equivalent to the flat case, so 
the minimum value of the surface frcc-cncrgy density in 
this region, ////, will be apfwi^ ~ 2i?c/i)/ cos a, with 
corrections of order of exp(— 77/^). Next Sections will be 
devoted to the evaluation of the minimum values of the 
surface free-energy functionals at the remaining regions, 
// and ///. 
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where Ci is the contour parallel to the boundary between 
regions / and /// between a wedge and apex (sec Fig. 

m- 

As argued in Refs. [l3l - [l5l [26j . the presence of cusps in 
the substrate induces the formation of disclination lines, 
which arc the responsible of the non-Berreman scaling of 
the clastic contribution to the surface free-energy density. 
In this Section, we will complete that analysis, evaluating 
the next-to-leading contribution. 



A. Singular contribution 

The solution 9{r) to the Laplace equation subject to 
the boundary conditions mentioned above can be split 
into two terms: a singular contribution 9s(v) due to the 
periodic array of disclination lines nucleated at the ridges 
of the substrate, and a non-singular contribution 0„s(r). 
A representation of the singular contribution for each tex- 
ture is given by [26j : 
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III. EVALUATION OF // 

The variations of 5* are restricted to the neighborhood 
of the substrate of a width typically of order of ^, and 
inside the defect cores. So, in region /, S takes the bulk 
value 5*;, everywhere [26j . Thus, the surface free-energy 
functional to minimize in region / reduces to a Frank- 
Oseen functional: 
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, q{y-Lsma) 
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where the integration is restricted to region /, AT is the re- 
duced elastic constant, and 9 is the orientation field. The 
minimization of Eq. ([7]) yields to the Laplace equation 
for 9, V^0(r) = 0. In the far field, we impose Dirichlct 
boundary conditions liiriy^oc ^(r) = aoo, where aoo = 
for the texture and aoo ~ 7r/2 for the A^" texture. 
Along the contours x = ±A/2, we should impose peri- 
odic boundary conditions. However, we impose instead 
Dirichlct boundary conditions 9 = aoo along these con- 
tours, as we know from the full LdG model minimization 
that these are the conditions satisfied by the mean-field 
solution [26j . Finally, we assume strong anchoring con- 
ditions along the boundary between regions / and ///: 
9{x, y = x^ tana/|a;|, z) = aoo + {x/\x\){a — aoo)- As dis- 
cussed above, this condition will be accurate if wL ^ 1. 
Their contribution to the free-energy density, //, comes 
from a contour integration of the mean-field solution via 
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Their contribution to the free-energy density, /|, comes 
from a contour integration of these solutions: 
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In Rcf. it was estimated the large-L behaviour of Eq. 
PT|) . leading to the non-Berreman term IC{a)q\n{L/ Rc) 
term. However, after some algebra it is possible to obtain 
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explicitcly /| from the solutions Eqs. ((9)) and (|T0)) via 
Eq.dni) as: 
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where /C(a) depends on the texture and substrate geom- 
etry as Eq. dH]) and we neglected terms of order q^{RcY 
and qrj. As the dependence on the nematic texture comes 
from /C(a), // will be minimum for the N-^ (resp. N^^) 
texture for a < 7r/4 (resp. a > n/A). 

Two remarks are pertinent at this point. First, we note 
that, although ff may depend on the explicit representa- 
tion of the singular solution, the leading non-Berreman 
contribution is independent on this representation. The 
reason for this fact is that this leading contribution arises 
from the behaviour close to the wedges and apexes of 6s, 
which must asymptotically approach to the correspond- 
ing to an isolated disclination line [2^. Secondly, the 
next-to-leading contribution gives a first contribution to 
B^ajw), which we will denote as Bj s(a), given by the 
expression: 



Bi,s{a) =IC{a) In 
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B. Non-singular contribution 



a tan a 



(13) 



The non-singular part of the mean- field solution, 6ns, 
is solution of the Laplace equation V^6ns = 0, subject 
to the boundary conditions 6ns = in the far field, i.e. 
y — >■ oo and along the boundaries x = ±A/2. On the 
other hand, 6ns = a — 6s along Ci. We do not have an 
explicit expression for 6ns (however, there is an implicit 
expression via a Schwarz-Christoffel transformation [l3l |. 
see below), so we have to resort to numerical methods. 
We have used two different techniques: a finite element 
method, analogous to the method outlined in the Sec- 
tion II to solve the LdG model, but minimizing instead 
the functional Eq. ([7]) subject to the boundary condi- 
tions for 6ns mentioned above; and as an alternative, the 
boundary element method [sj, [s^l • In this technique the 
solution 6ns inside the region I can be written as: 

6ns{r) = £^ds|^[z/(s) • V,0„,(s)]G(s,r) 

-6ns{s)Hs)-V,G{s,r)]^ (14) 

where the contour integral over the boundary dl of region 
/ is counter-clockwise, i/{s) is the outwards normal to the 
boundary at s and G{s, r) is the fundamental solution of 
the Laplace equation G(s, r) = — In (|s — r|) /27r. As we 
impose Dirichlct boundary conditions, the second term 
in the right-hand side of Eq. (ITil) is known. On the 



other hand, the normal derivative of 6ns at the boundary 
is obtained by solving the integral equation [s^ [13] : 

/ ds[i.(s).VA,(s)]G(s,r)==^^ 
Jai ^ 

+ (f ds0„,,(s)[z/(s). VsG(s,r)] 

Jdl 



(15) 



where now r e dl. In order to solve Eq. (fTS)) . we dis- 
cretize the boundary in a set of straight segments (the 
boundary elements). We use the constant boundary ele- 
ment approach [l^], so we assume that both 6ns and its 
normal derivative are constants along each boundary el- 
ement. Introducing this approximation to Eq. (|15|) . we 
obtain a set of linear algebraic equations for the normal 
derivatives of 6ns- Once we solve this set of equations, 
and introducing the same approximation in Eq. (1141) we 
obtain the non-singular orientational field 6ns inside re- 
gion /. 

Once evaluated 6ns , its contribution to the surface free- 
energy density can be obtained from: 
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where the first result is more appropriated for the bound- 
ary element technique, while the second is more appro- 
priate for the finite element method (note that the last 
term in the second result can be evaluated numerically 
with high accuracy by standard methods as we know an- 
alytically 6's(r)). 

The numerical minimization is performed in the cell 
shown in Fig. [2l In order to minimize finite-size effects, 
the cell height B is taken to be at least four times the 
value of L (note that 6s decays exponentially to oioo for 
y 2Lcosa/n). We checked that the value of A/"^ 
is independent of A, as expected, so the non-singular 
contribution to B{a,w), Bi^ns{oi), is related to via 
Bi^ns{oi) = A/I"*. The numerical results are shown for 
both the and A^H textures in Fig. [31 The agree- 
ment between the results obtained from the finite element 
method and boundary element method is excellent, as for 
the evaluation of 6'„s(r), which takes non- negligible val- 
ues only above the contour Ci and vanishes close to the 
wedges and apexes (see Fig. |4]). Our results show that 
the non-singular contribution corresponding to the A^^ 
(resp. A^") texture is smaller than the contribution asso- 
ciated to the A^ll (resp. A^-'-) texture for a < n/i (resp. 
a > 7r/4). 



C. Exact evaluation of // 

In the previous paragraphs we have outlined how to 
obtain 6{r) — 6s{r) + 6'„s(r), and from that, to obtain 
// = // + as weU as Bj{a) = Bj,s{a) + B/,„,(a). 
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However, it is possible to obtain // directly without 
knowing the explicit form of 9(r) by using a Schwarz- 
Christoffel transformation [13, y3 : 



/ 



c 



(C + l)l/2-a/'r^2a/77(^ _ iy/2~a/7r 



a I a 3 



where z = x + iy, 2Fi{a,b,c, z) is the Gauss hyperge- 
ometric function and C, C" and C" complex constants. 
This conformal transformation maps the minimization 
cell for — >■ oo into the upper half C^-plane (see Fig. 
[2]), transforming the origin into the origin and the edges 
z = ±Lcosa + zLsina into C = ±1, respectively. These 
conditions fix the values of C and C", so the Schwarz- 
Christoffcl transformation finally reads: 



£i r 

TT / 

a 1 
^'2 



TT / 

a 3 
^' 2 



(19) 



where r(a;) is the gamma function. Eq. (|19p can be 
formally inverted, so = (^{z/L;a) = x' {x/ L,y / L;a) + 
iy'{x/L, y/L; a). We will consider the limit 77/L — 0, and 
small but finite. In this approach, the boundary 
of zone // becomes under the Schwarz-Christoffel trans- 
formation the real axis in the C-plane, rounded around 
C = and ±1. The expansion of Eq. ([T9l) around these 
values show that the circles of radii Rc around the origin 
and the edges z = ±L cos a + iL sin a map into circles (up 
to corrections of order of {Rc/ L)^) of radii ei for C = ±1 
and £2 for C = given by: 
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of X and y, so we cannot give an analytic expression for 
6{x, y). However, we can evaluate exactly // since: 



As the Schwarz-Christoffel transformation is conformal, 
and 9{y) is harmonic, we are going to find the solution 
to the Laplace equation in the half ^-plane, 9, subject 
to the boundary conditions 9 = Uoo for \x'\ > 1 -I- ei, 
and 9 = aoo + {x' /\x'\){a — aoo) for €2 < \x'\ < 1 — ei. 
The solution 9{x' , y') in the image of the region / on the 
C-plane, /', is given by [isj : 



(a;', y') = ctoo H — arctan 
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From the solution 9, we can obtain 9{x, y) = 
9[x'{x/L,y/L\a),y'{x/L,y/L;a)]. We note that we do 
not have an explicit expression for x' and y' as functions 
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where B is the image in the ^-plane of the boundary of 
the zone /. After an straightforward calculation, we get 
the expressions for // and Bj{a): 
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Fig. [S] shows an excellent agreement between the theo- 
retical prediction for Bj{a) from Eq. ([M)) and the re- 
sults obtained in the previous subsections, except close 
to a = tt/2. The latter may be due to numerical uncer- 
tainties in the evaluation of Bj^ns, since either it diverges 
{N-^ texture) or vanishes (A^" texture) in that limit. As 
happened for Bj s, the dependence on the texture comes 
from IC{a), leading to the same conclusions about the 
relative stability of the nematic textures with the angle 
a, in agreement with the results obtained in Ref. jisj . 



(20) 



IV. EVALUATION OF /// 

In Ref. [l^ it was observed that there is a dependence 
on w of the next-to-leading contribution to the surface 
free-energy density. However, in the previous Section we 
have shown that the contribution from region I only de- 
pends on a. So, we anticipate that this dependence comes 
from region II, where inliomogoneities of both the ne- 
matic order parameter S and orientational 9 fields are 
observed. 

The free energy of region // can be evaluated as the 
sum of the contributions of the regions around each cusp 
(either wedge or apex). In each of these regions and if 
Rc <C L, we anticipate that the orientational field far 
from the cusp behaves asymptotically as that of the iso- 
lated disclination line which has been nucleated at the 
substrate wedge or apex: 9 ~ /(/), where / is the topo- 
logical charge of the disclination line and (r, (p) are the 
polar coordinates taking as origin the cusp. The val- 
ues of / are fixed by the (strong) anchoring conditions 
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on the substrate and depend on the nematic texture: 
Ii = — a/(7r/2 — a) (resp. Ii = +1) at the wedge bottom 
and I2 = a/(7r/2 + a) (resp. ~{tt/2 - aV(7r/2 + a)) at 
the apexes for the N-^ (resp. iVil) [ll,[23. However, S 
decreases from the bulk value far from the origin and van- 
ishes as r — > in order to remove the free-energy singu- 
larity associated to the defect core. Although we cannot 
solve analytically this problem, we can estimate using an 
ansatz the free energy associated to this combined distor- 
tion of 9 and 5, as shown in the Appendix. However, we 
resort here to a full minimization of the LdG model Eq. 
([T]) restricted to region //. At r = i?c, we impose Dirich- 
let boundary conditions to Q, where S takes the value 
corresponding to the order-parameter profile correspond- 
ing to a planar wall at the distance between the boundary 
point and the closest substrate, the biaxiality parameter 
B = and 6 = a^o + I{4> ~ ""/S)- Alternatively we used 
free boundary conditions v ■ VQ = 0, leading to similar 
results. Fig. [5] shows typical textures obtained after min- 
imization. At distances r larger than a few correlation 
lengths, S decays to the bulk value except in the neigh- 
bourhood of the substrates, where takes approximately 
the value corresponding to the flat substrate profile. On 
the other hand, the orientational field deforms contin- 
uously in order to satisfy anchoring conditions on each 
side of the wedge or the apex. So, we anticipate that the 
contribution to the surface free-energy density /// will 
scale with Rr as: 



Mil = 4cr„^,i?c + /C(a) In + Bn 



(25) 



The first contribution arises from the inliomogeneities of 
S close to the substrates, and the second one from the 
asymptotic behaviour of for r > 1: by using the Frank- 
Oseen functional Eq. (O, and taking into account that 
V^? ~ lu^/r (where is the azimuthal unit vector), 
we find that the free-energy contribution is proportional 
to {K/2)I^ A(j) In Rc, where Acp is the opening angle of 
the wedge or apex. From this result we get the second 
contribution in Eq. ([^5)) by noting that /C(a) = Klf{TT~' 
2a)/2 + Kl2{n + 2a)/2. The remaining contribution will 
give the core free energy per unit length associated to the 
disclination lines. 

Expression (|25|) can be used to extract Bjj from the 
numerical minimization. So, fixing the value of a, and 
for each w (between and 1.5) we considered a range of 
values of Rc between 10 and 90. The minimization was 
performed by the mesh-adaptive finite-element method 
used previously [32|, taking t = 1 and k = 2 (this choice 
is motivated to compare with results reported in the lit- 
erature [2^, H^). After substrating AamoRc to Xfii, the 
numerical results clearly show a logarithmic dependence 
on Rc, with a slope approximately equal to ]C{a) (see 
left panel in Fig. [7]). Next step was to substract the sub- 
dominant contribution /C(a)lni?oo- Now the numerical 
results are nearly independent of Rc (see right panel in 
Fig. [T]). The value of Bu is estimated as the mean value 
of these results, with an errorbar given by the dispersion 
of the numerical data around the average. 



V. DISCUSSION 

The results obtained in the previous Sections can be 
combined as follows: 
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up to corrections of order of q^. Consequently: 



(27) 



B{a, w) ~ Bj{a) + Bji{a, w) = 

1 a \ . / + a 
2^ n 



-IC{a) 



In 



-"HI- 

+Bii{a, w) 



(28) 



This is the main result of our paper. Now we can check 
this prediction by comparing these results with those ob- 
tained from the full minimization of the LdG model of a 
nematic in contact with a sawtoothed substrate [2^. Fig. 
|8] shows the comparison between the results reported in 
the Ref. [IBl (r = 1, k = 2) and the calculated ones 
in this paper for a = tt/6 {N-^ texture) and a = tt/3 
(A^ll texture). For large w, the agreement is good, al- 
though our results slightly overestimate those from the 
full minimization. On the other hand, for small w the 
curves obtained from the full minimisation converge to- 
wards our prediction as L increases. So, our approxima- 
tion is accurate even for moderate values of L, despite 
the assumptions involved in our approach. Consequently, 
the scheme considered in this paper is an alternative to 
the full minimization technique, which is quite expensive 
from a computational point of view, when wL > 1. 

Finally, it is interesting to note that the elastic contri- 
bution to the surface free-energy density, Eq. ((27| . can 
be expressed as follows: 



A/, 



elastic 



2 X 



2.i,'^l^lnlM 
2 2 L 



+ Bii{w,a) (29) 



where we note that IC{a) ^ -4ttK{Ii/2){I2/2), with Ii 
and I2 being the topological charges associated to the 
disclination lines at the wedges and apexes, respectively. 
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The characteristic length 7(a), which absorbs the contri- 
bution to B from region /, is defined as: 



7(a) = r 
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(30) 



This lengthscale decays continuously from tt/2 for a — > 
to 1 for a —s- tt/2. We can understand the first contribu- 
tion in Eq. ((29)) as the interaction between the disclina- 
tion line at x = with half its topological charge and the 
apexes at a; = ±Lcosa, again with half their topologi- 
cal charges [s^ . As the absolute values of the topological 
charges are always smaller or equal than 1 , the nucleation 
of disclinations at the cusps of the surface is favourable 
with respect to bulk disclinations (which only may have 
half-integer values). With this intepretation, 7(a) can be 
used to define an effective core radius 7(a), and the sum 
of the core energies is given by Bii{a,w). 



VI. CONCLUSIONS 

In this paper we have analyzed elastic contribution to 
the surface free-energy density for a nematic in contact 
to sawtoothed substrates in the strong anchoring regime, 
i.e. wL > 1. We have extended the analysis done in Rcf. 
[26| . So, in addition to the leading contribution propor- 
tional to — qlnq, with q being the substrate periodicity 
wavenumber q, we have characterized the next-to-leading 
term. This term has two contributions: one associated to 
the deviation of the orientational field with respect to the 
contribution of the array of disclination lines nucleated 
at the cusps of the substrate, and the core free-energy 
associated to them. We anticipate that our analysis can 
be generalized for other substrate shapes when the ne- 
matic texture presents topological defects induced by the 
structure, as, for example, in crenellated substrates [ssj . 
Furthermore, our results may be used to predict accu- 
rately the location of first-order wetting transitions in 
nematic liquid crystals in contact to general substrates 
[28I [2^, HI]. However, we must note that our analysis is 
restricted to bulk nematics liquid crystals with in-plane 
deformations in presence of grooved substrates. The ef- 
fect of the substrate structure in the nematic texture in 
partially filled configurations (i.e. with an isotropic fluid 
in bulk) , the influence of twist nematic deformations [l^ 
or the effect of the structure in the z direction deserve 
further study, but this is beyond the scope of this work. 



Acknowledgments 

The authors wish to thank Dr. P. Patricio and Dr. 
N. M. Silvestre for very stimulating conversations and 
technical advice on the finite-element method, and Prof. 
M. M. Telo da Gama, Prof. L. F. RuU and Prof. A. O. 
Parry for enlightening discussions. We acknowledge the 



support from MICINN (Spain) through Grant FIS2009- 
09326, and Junta de Andaluci'a (Spain) through Grant 
No. P09-FQM-4938, both co-funded by EU FEDER. 



Appendix A: Ansatz for the evaluation of /jj 

In order to estimate the free energy associated to the 
region //, we can make use of an ansatz. As argued in the 
paper, the nematic order parameter S must vanish at the 
cusps of the substrate to avoid the divergence of the free 
energy, but it must converge to the bulk value far from 
the wedges and apexes. So, in the region around each 
cusp, we can suppose that S depends only on the radial 
distance r from the wedge or apex. In particular, we use 
the following ansatz [s^]: S{r) = Sb{l — e~^), where St 
is the bulk nematic order parameter and /3 is a length- 
scale to be determined later. We will assume that the 
biaxiality parameter B = everywhere. Regarding the 
orientational field 9, we assume that has the expression 
corresponding to a disclination line placed at the origin: 
9 = ctoa + 1(4' ~ ■"'Z^), where / is the winding number 
associated to the disclination line, and aoo is the far-field 
value of the orientational field (0 for the 7V^ texture, 
and tt/2 for the A'^" texture). Substituting this ansatz in 
the LdG functional Eq. ([1]), this expression reduces to a 
frmction of /3: 



A/ « ^(tt T 2a)f3' - 2w{Rc ^ /?) + y ] 9/2(7r T 2a) 



+ ln(^ 

V2/3 



16^ ^ 8 V (/-I) 



=Fl8/sin(2a) - 3 [ sin(-2/(7r T «) + {I - l)7r) 
-sin(-2/(±a) + (/-l)7r)])| (Al) 

where we have considered the case r = 1 and k ~ 2 (for 
other values, we can get straightforwardly analogous ex- 
pressions). In this expression, the upper sign corresponds 
to the wedge situation (so / = /i), and the lower one to 
the apex (so I ~ I2); and 7 is the Euler constant. Mini- 
mizing Eq. (|A1|) with respect to /3, we obtain the optimal 
value for this lengthscale: 



/3 



72 



w 



13 (tt =f 2a) 




13 (tt t 2a) 



648 



P (A2) 



Substituting Eq. (|X2|) into Eq. (PT|) . we get the free 
energy estimate associated to each region around a cusp, 
and provides an upper limit to the real value of A/. 

In order to have an estimate of -B//, we substract to the 
obtained values of A/ through this ansatz the surface and 
elastic contributions —2wRc and (A'/2)/^(7r=F2a) InRc = 
(9/7)/^(7r=F 2a)lni?c, respectively. After this, the esti- 
mate of Bjj is obtained as the sum of the results obtained 
for the wedge and apex. Fig. [S] shows the comparison 
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of this estimates with the results obtained from the full although our ansatz overestimate the value of Bjj, and 
minimization, for a = 7r/6 [N^ texture) and a = 7r/3 qualitatively the dependence of i?// with w is recovered. 
(TVil) texture. They have the same order of magnitude, 
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FIG. 1: Schematic picture of the geometry of the system, 
characterized by the side length L and the angle a. The 
different regions I, II and III are outlined. See text for 
explanation. 



H 




a) b) 

FIG. 2: Left panel: Minimization cell for the evaluation of 
fi,ns- Right panel: Mapping of the minimization cell under 
the Schwarz-ChristofFel transformation in the limit H — oo. 
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a(rad) 

FIG. 3; Plot of \fi^ns/K as a function of a for the A'"^ and 
A^" textures, by using the finite-element method (F.E.M.) and 
the boundary element method (B. E. M.). 




FIG. 4: (Color online) Plot of 9ns as a function of x and 
y obtained from the finite-element method for a = n/Q and 
L = 0.25. For the sake of clarity it is only represented in half 
a minimization cell < a; < L cos a. 
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FIG. 5: Plot of Bi/K.[a) as a function of oi. The continu- 
ous broad line corresponds to the exact expression, and the 
symbols correspond to the estimates from the numerical data 
for fi,ns (from the boundary element method): circles for the 
A''^ texture and diamonds for the A''" texture. Inset: plot of 
— Bi/K as a function of a. The meaning of the symbols is 
the same as in the main panel. 



V' 

FIG. 6: (Color online) Typical textures obtained from the 
full minimization of the LdG function in region // for r = 1, 
K, = 2, w — 1.0, 1/ = 90 and a = 7r/6 (upper panels) and 
Of = 7r/3 (lower panels). The color map corresponds to the 
reduced nematic order parameter field S — S/St, and the 
segments correspond to the nematic director field n. 
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FIG. 7: Left panel: representation of Xfu — AgnwRc as a 
function of Rc for w = 1 and: a — n/G (squares) and a — n/S 
(circles). The straight lines correspond to the logarithmic 
regressions of the numerical data. Right panel: representation 
of Bii{a, w) as a function of Rc. The meaning of the symbols 
is the same as in the left panel. 




FIG. 8: (Color online) Comparison between the values of 
B(a, w) obtained from the full minimization of the LdG model 
[2b | and the results obtained in the present work. Symbols 
correspond to B{a, w) obtained from the full minimization of 
the LdG model for different values of L and q = 7r/6 in the 
N''- texture (black symbols), and a = tt /3 in the A/'" texture 
(orange or light gray symbols). The wide blue (black) and 
wide red (dark gray) lines with errorbars correspond to the 
predictions from this work for a — n/6 and n/3, respectively. 




FIG. 9: Comparison between the values of Bii{a,w) ob- 
tained from the full minimization of the LdG model and the 
ansatz (see text). Symbols correspond to Bii{ce,w) obtained 
from the full minimization of the LdG model for a = 7r/6 in 
the N'^ texture (circles), and a = 7r/3 in the A'^" texture (up 
triangles). The wide continuous and dashed lines correspond 
to the estimates from the ansatz for a — 7r/6 and a = vr/S, 
respectively. 



